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I We introduce and study a family of quantum fields, associated to 5-interactions in one di- 

■ mension. These fields are analogous to holonomic quantum fields of M. Sato, T. Miwa and 
M. Jimbo. Corresponding field operators belong to an infinite-dimensional representation of 

I the group SL(2,M.) in the Fock space of ordinary harmonic oscillator. We compute form fac- 

1/^ ' tors of such fields and their correlation functions, which are related to the determinants of 

■ Schroedinger operators with a finite number of point interactions. It is also shown that these 
I determinants coincide with tau functions, obtained through the trivialization of the det*-bundle 

T— I ' over a Grassmannian associated to a family of Schroedinger operators. 

vn : 
o ■ 

1 Introduction 

'-^ I The study of holonomic quantum fields ^I] has led to important advances in both integrable quantum 
field theory and analytic theory of linear differential equations. In physical language, these fields 
represent a particular case of the canonical Bogolyubov transformations (see, for instance, 0), which 
makes possible the exact computation of their form factors and correlation functions. 
^ , The main physical examples of holonomic quantum fields are given by the exponential fields of the 
sine-Gordon theory at the free-fermion point (SG//) [SIEIII21 and order /disorder variables of the two- 
dimensional Ising model [HI E] . Both models possess an underlying free-fermion structure. However, 
their observables are nonlinear in terms of free fields (they are represented by ordered exponentials 
of fermion bilinears) and have nontrivial braiding relations with them. This nonlinearity leads to 
interesting correlation functions. 

Correlators of holonomic quantum fields are usually called tau functions. In the special case of the 
SG//-theory they have the meaning of determinants of Dirac operators with branching points on the 
Euclidean plane. An attempt to give a geometric definition of such tau functions was made in |10j . 
It was based on the approach, developed in the article [H], where the tau function of the Schlesinger 
system was related to the determinant of a singular Cauchy-Riemann operator. The idea of PU] was 
to consider a family A of Dirac operators, parametrized by the coordinates of branchpoints, and to 
associate to each of these operators a subspace of boundary values of local solutions of the Dirac 
equation. These subspaces are then embedded into an infinite-dimensional grassmannian. One can 
construct d la Segal- Wilson JH] the det*-bundle over this grassmannian and its canonical section a. 
Next, using the Green function of the Dirac operator, one may endow the det*-bundle with another, 
trivializing section. This latter allows to identify a with a (tau) function on A. 
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The main goal of the present paper is to explain the concept of holonomic quantum fields and 
the above definition of the tau function with a simple example. It appears that these two notions 
naturally emerge in the calculation of the resolvent of the Schroedinger operator with 5-interactions in 
one dimension. Though such operators and their resolvents have already been extensively discussed 
in physical and mathematical literature (see, for example, the monograph pp and references therein), 
holonomic quantum fields allow to examine this fairly classical subject from a new point of view. 

Moreover, we believe that the ideology developed in this paper applies as well to certain unsolved 
quantum mechanical problems, the most appealing one being the computation of the fermionic 
vacuum quantum numbers, induced by a finite number of magnetic vortices in 2 + 1 dimensions. In 
practice, the latter problem reduces to the calculation of the resolvent of the Dirac hamiltonian with 
point sources. At present, the answer is known only in the case of a single vortex |14j. It seems, 
however, that the multivortex resolvent can be obtained from the correlation functions and form 
factors of certain Bogolyubov transformations, generalizing exponential fields of the SG//-theory. 

This paper is organized as follows. After introducing basic notations and terminology in Sec- 
tion 2, we turn in the next section to the calculation of the resolvent of the Schroedinger operator 
with 5- interactions. It is expressed (by the formula (j3.3|) ) through the ratio of correlation functions 
of certain local fields in the ID quantum field theory of free massive real bosons. These correlation 
functions are computed in the lagrangian approach by an auxiliary integration method. They com- 
prise the fields of two types: the free ones and interacting fields, associated to delta- sources. These 
latter represent the simplest prototypes of holonomic quantum fields, since the operators, correspond- 
ing to them in the hamiltonian picture, realize certain canonical Bogolyubov transformations of the 
Heisenberg algebra. In Section 4, we compute form factors and correlation functions of the fields, 
which correspond to more general Bogolyubov transformations and belong to an infinite-dimensional 
representation of the group SL{2, M) in the Fock space of harmonic oscillator. Section 5 is devoted 
to the definition and calculation of the tau function. It is obtained through the trivialization of the 
det*-bundle over a finite-dimensional grassmannian of boundary conditions, associated to a family 
of Schroedinger operators with point interactions. Appendix A contains several formulas of real and 
complex gaussian integration. We conclude with a brief discussion of possible generalizations and 
open questions. 

Acknowledgements. The author is grateful to A. I. Bugrij, M. Z. lorgov, V. N. Roubtsov, 
V. M. Shadura and Yu. O. Sitenko for stimulating discussions and useful comments. 



2 Schroedinger operator without point interactions 

Let us first consider Schroedinger operator without point interactions, 




acting on functions from the Sobolev space if^(M) as the second derivative. Its resolvent {L — E) ^ 
is given by the integral operator with the kernel 

-m\x-y\ 

which can be easily evaluated by Fourier transformation. Here we have introduced the notation 
E = —m? and assumed that Rem > 0. 

Note that for real < the resolvent kernel Ge{x, y) coincides by construction with the two-point 
correlation function in the relativistic euclidean quantum field theory of free massive real bosons in 
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one dimension. In particular, if we define the action 



So[<p] = \ j dx ip{x) (^-^ + rri^ ip{x), (2.2) 

— oo 

then GE{x,y) can be formally written through the ratio of two functional integrals: 

Ge{x, y) = {ip{x)ip{y)) = j ^s^M • (2-3) 

The theory, described by the action ()2.2p . admits two natural interpretations: 

• If we interpret our single dimension as space, then the action ()2.2|) coincides with the energy 
functional of the infinite string in a parabolic well. The integral in the denominator of (j2.3j) 
represents string partition function, and two-point correlator {ip{x)ip{y)) is the thermodynamic 
average of the product of transverse coordinates of two different points of the string. 

• On the other hand, the action (j2.2j) describes the dynamics of harmonic oscillator in imagi- 
nary time. In this setting, correlation function {(f{x)(p{y)) may be interpreted as the vacuum 
expectation value of the ordered product of coordinate operators at two different times. 

In order to introduce several important notations, let us briefly recall the hamiltonian approach to 
oscillator dynamics. Fields ip{x) and 7r(x) = 7^^;^ here become operators, obeying the commutation 
relation [0, tt] = i. The hamiltonian, being expressed in terms of (f and vr, has the form 

where the constant term is subtracted for future convenience. Imaginary time evolution of an arbi- 
trary operator O is given by the equation 

d{x) = e-^"C'(0) e^". (2.4) 

It is customary to define the creation-annihilation operators 

m ( ^/^^ , Z7r(0)^ I /"^/-/n\ ^7r(0) 



^(0) + ^^ , aT = ./-U(0) 

satisfying canonical commutation relation [a, a^] = 1. The hamiltonian can then be rewritten in 
terms of these operators as H = ma^a. Vacuum vector |0) is fixed by conditions a|0) = 0, (0|0) = 1. 
The operators that we wish to consider act in the Fock space JF, spanned by the orthonormal vectors 

\k) = ^\0), fc = 0,l,2..., (2.5) 



constituting the set of hamiltonian eigenstates: H\k) = km \ k). 

The computation of correlation functions of local fields in the hamiltonian approach is equivalent 
to the calculation of form factors, i. e. matrix elements of the corresponding field operators in the 
orthonormal basis of eigenstates of H. For instance, form factor expansion of the two-point correlator 
((9i(xi) 02(3^2)) is written as 

(Oi(xi)02(x2)) = (0|di(a;i)a2(x2)|0) = (0|ai(0) 6-^(^^-^^)02(0)10) = 
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= J2 (0| C'i(O)lfc) (^102(0)10) e-'^^^-^-xi)^ ^2.6) 

fc=0 

where we have assumed that X2 > xi. One can obtain analogous expressions for the multipoint 
correlation functions, simply using the formula ()2.4|) . and inserting the appropriate number of times 

oo 

the resolution of the identity operator, 1 = ^ \k) {k\, into the correlators. 

Example. It is easiest to see how this scheme works on the example of the two-point correlation 
function {ip{x)ip{y)). Note that since the operators (p and vr are given by 




the only non-zero form factors are 



{k + l\0{O)\k) = {k\0{O)\k + l) = 



+ l|7r(0)|A;) = -(A;| 7r(0)|A; + 1) = z 

that is, the numbers of in- and o'ut-particles should differ by 1. Therefore, in the expansion over 
intermediate states ()2.6|) for {if{x)(f{y)) only the terms with k = 1 will remain, reproducing thus the 
formula (EH). 





3 Introducing (^-interactions 

Let us now see what happens if instead of L we consider Schroedinger operator with a finite number 
of 5-interactions, 

Lay = + V{x), Vix) = J2^^ - • 

1=1 

The most known way of calculating the resolvent {Lay — E)^^ is to expand it formally in a se- 
ries in powers of V. Summing up this series, one obtains a compact expression for the resolvent 
kernel Ge,v{Xj u)'- 

N 

GE,v{x,y) = GE{x,y) - ^ GE{x,ai) Ufj^ GE{aj,x), (3.1) 

Here GE{x,y) denotes the unperturbed resolvent ()2.H) and the matrix U is defined as 

Uij = ^ 6ij + GE{ai, aj), i,j = l,...,N. 

An alternative simple proof of this result follows from field-theoretic considerations. Again, for 
real negative E = —m?' the resolvent GE,v{x,y) coincides with the two-point correlator in the one- 
dimensional quantum field theory, described by the action 

oo 

Svlv] = \ j dx ^{x) (^-^ +m^ + V{x)^ ^{x) = SM + SM. (3.2) 

— oo 

4 



where 

N 



Thus we have 
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1=1 



G's,y(x,i/) = {^{x)ip{y))v = jp^ ^-sy[^] • 

The action ()3.2|) reproduces the energy of a one- dimensional string with N masses attached to the 
points Oi, . . . , a^r. We will assume for definiteness that these masses are all positive. 

The idea is to represent the factor e~^'^^''^\ appearing in the functional integrals, as a gaussian 
integral over N auxiliary variables fj,i, . . . , fi^: 



(27r) 



N/2 



— — / ... / c//ii . . . rf/i^ exp <^ - ^ /i//V;- + ^ ¥?(aj) \ 

lli=iV^ijL jL I i=i i=i J 



Now one can make use of the formula ()A.2|1 of the Appendix. The matrix A in our case has continuous 
indices, and its inverse should be replaced by the resolvent {L+m?)~^ . The matrix B is diagonal, 

Bki = ^ Ski, k,l = 1, . . . ,N, 

and the matrix C has one continuous and one discrete index: 

Cx,k = —i S{x — Ofc), k = 1, . . . , N, X G M. 

Substituting these matrices into the relation ()A.2|) . we immediately obtain the formula ()3.1|) for the 
resolvent kernel GE,v{x,y). 



Remark. Schroedinger operators with 5-interactions at the points ai, . . . , oat are rigorously defined 
as the eleme 
the domain 



as the elements of an A^-parametric family of self-adjoint extensions of the operator La = with 



domL, = {/ G E\^) I /(a,) = 0, j = 1, . . . , AT} . 

This operator has deficiency indices (A^, A^) and thus the whole family of its self-adjoint extensions 
depends on A^^ parameters that can be incorporated into a hermitian N x N matrix B. Similarly to 
the above, one can associate to a general self-adjoint extension La^B the action 

°? AT 

Sb[^,0\ = / dx (p{x){L + m'^)(p{x) + ^ (^(ai)5.jV(«i)- 
The kernel Ge,b{x, y) of the resolvent {La^B + fn'^)~^ coincides with the two-point correlation function 



J Dip Dip e'^B[v>,^] 



It can be easily evaluated, using the above method of auxiliary fields and applying the formulas of 
complex gaussian integration from the Appendix. The result has the form (j3.1|) : the only difference 
is that the matrix U now depends on B: 

Uij = Bij + GE{cLi, aj). 
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In what follows, however, we will be dealing only with local point interactions. 

Let us consider yet another approach to the calculation of the resolvent GE,v{x,y). On the one 
hand, it is equal to pair correlation function in the theory with the action ()3.2j) . However, it may 
also be expressed through the ratio of certain multipoint correlators in the theory, described by the 
unperturbed action ()2.2|) . Namely, one can write 

Ge,v{x, y) = {ip{x)ip{y))v = , . -z^, — rr , (3.3) 

where the local fields Oy(a) are defined as 

Oy{a)=e^v\-\vv\a)Y (3.4) 

It should be pointed out that the correlation function (Cv'j(ai) . . . Cy^(aAr)), standing in the 
denominator of (|3.3p . is equal to the ratio of partition functions of the string with and without 
attached point masses. It can be formally expressed through the determinants of Schroedinger 
operators: 

(O..(a0...a.M) = ^/Jg^. (3.5) 

Repeating the trick with auxiliary fields and using the formula ^A.l\ of the Appendix, one may also 
compute this correlator exphcitly: 

/ \ -1/2 

(Oy,(ai)...Ow(a^)) = f det 5^^ + ^Wj G E{ai,a,) \ , ^, j = 1, . . . , AT. (3.6) 

Consider now the fields Oy in the hamiltonian picture. The operator Oy = exp | — ^ V(p^^ has 
remarkable equal-time commutation relations with the operators and if: 

dy{o)m - movio) = 0, (3.7) 

dv(0)7r(0) -7r(0)C>y(0) = -iV0{O)dyiO). (3.8) 

This is a manifestation of the fact that any correlator f{x) = {. . . Oy{a)(f{x) . . .), being considered 
as a function of x, is a local solution of the Schroedinger equation (L — £")/ = 0, satisfying at the 
point a the boundary condition 

/(a + 0) - /(a - 0) = 0, f'{a + 0) - /'(a - 0) = Vf{a). 

One can also rewrite the commutation relations ()3.7p - ()3.8p in terms of the creation-annihilation 
operators: 

o.,o) ( ) ^ ( i-.)(:0- f^.^) 

Therefore, Oy realizes a Bogolyubov transformation, i. e. linear transformation of the Heisenberg 
algebra, preserving canonical commutation relation [a, a^] = 1. The formula ()3.9j) determines the 
operator Oy almost completely. More precisely, Oy is fixed by ()3.9j) up to a constant numerical 
factor. 

In the next section, we will compute form factors and correlation functions of the fields, correspon- 
ding to more general Bogolyubov transformations. It should be emphasized that we will use only 
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the relations of type (j3.9p and no reference will be made to the explicit formula (|3.4p . The reason for 
doing so is that it is not always clear, which operator should be associated to a given field and vice 
versa. This is the case, for instance, in the analysis of ^'-interactions. Even more severe difficulties 
arise in some two-dimensional problems: the fields, realizing relevant Bogolyubov transformations, 
being themselves local, are not mutually local with the free fields. The most known example of this 
kind is given by the exponential fields in the SG//-theory. 

4 Form factors and correlation functions 
of Bogolyubov fields 

Consider a linear transformation of the creation-annihilation operators 

a \ ( ^ \ _ ( ^ P \ 1 



It preserves canonical commutation relation iff the real parameters a, f3, 7, 6 satisfy the condition 
a6 — jS'-f = 1. In the following, this condition is assumed to hold. We want to represent A as a 
similarity transformation. Namely, we are looking for the invertible operator 0\, such that 



These operators are called Bogolyubov transformations, and the corresponding local fields will be 
called Bogolyubov fields in the rest of this paper^. The operators Oa are determined by ()4.2|) up to 
a constant, realizing thus an infinite-dimensional projective representation of the group SL{2, M) in 
the Fock space of our harmonic oscillator. 

In order to construct 0\ in terms of a and a\ one should examine the properties of basic 
Bogolyubov transformations: 



00 



n=0 



It is easy to check that induced linear transformations have the form 



General Bogolyubov transformation ()4.2j) is then given by 



- - - f 1 2 1 

Ox,^,u = PxQf^Ru = : exp j - A (a^) + /i a^a + - 



ua'}: (4.6) 



^Actually, there already exist several names for some special fields of this type, to mention only 'holonomic quantum 
fields', 'exponential fields' and 'monodromy fields'. All these names, however, do not reflect common structure of such 
fields, and are related to the particularities of different two-dimensional problems. That is why we believe that the 
name 'Bogolyubov fields' would be more appropriate. 
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Parameters a, /3, 7, 5 of the corresponding linear map A = Aji^Aq^Ap^ are related to A, /x and u by 
the following formulas: 

« = T^, ^ = -T^' ^=T^' S = l + ^I--^. (4.7) 

l + /i 1+yU 1+ fl 

Note that writing 0\ in the form ()4.fj|l . we adopt the convention {0\) (OIOaIO) = 1. This normal- 
ization will be used hereafter. Together with the relation ()4.2|) . it completely fixes the operator 0\. 

Remark. One-point function of the Bogolyubov field Oy, associated to a 5-interaction, is not equal 
to 1. It may be determined from the formula ()3.6p by setting = 1: 

(O.) = (1 + ^) . (4.8) 

It can also be computed in the hamiltonian approach, being rewritten as the vacuum expectation 
value 



{Oy) = {0\dy\0) = (0|e-^^^>) = Y.^(-^) (« + «^)'"|0)- 



(4.9) 



n=0 

, 2n 



Using Wick's theorem, one obtains (0| (a + a"!") |0) = {2n — 1)!!. Now remark that the Taylor 
expansion of the function 



°° f2n- 1)!! / X 



n=l 



l + = l + (4.10) 



coincides with (14. 9p . if we set a; = Thus we recover the formula (j4.8p . 

Parameters A, fi and u, which correspond to the field Oy, are determined from the comparison 
of (jSini) and (jOI). The result is 

V/2m 

A = /i = z/ = — . 4.11 

^ 1 + V2m ^ ' 



This leads to the following representation of the operator O 



V- 



We now turn to the calculation of form factors {k\0\^^^v\l) of the operator ()4.6|) . It is clear 
that such form factors will be non-zero only if the numbers of particles in in- and out-stsXe are 
simultaneously even or odd. It is convenient to use instead of {k\0\^^^v\l) the auxiliary variables 



^\VT\ (fclOA,^,.!/) = {Q\a^PxQ^Ru (a^)' |0). 



Using the relations ()4.3|) - ()4.5p in the last expression, one can pull the operator P\ through , and 
the operators Q^, through (a^) : 

Fk,i = (0|A(a + Aat)' Q^R, (a^)' |0) = (0|A(a + Aat)' [u a + a))' R,\Q) = 



Next, since we have 

(0|Pa = {Q\Q>. = (0|, QJO) = ^,|0) = |0), 

the variable F^^i may be rewritten as the vacuum expectation value of the product of certain linear 
combinations of the creation-annihilation operators: 

Fk,i = (0| {a + Xa^Y ('^(l + yu)~'a + (1 + yu) a^)'|0)- 

Wick's theorem allows to express this vacuum expectation value through the sum over all possible 
pairings between the linear combinations. There will be only three types of such pairings: 



(0| (a + Aa^) |0) = A, (4.13) 

{0\(^u{l+ fi)-^a + {l + fi)a^y\0) = u, (4.14) 

(0| (a + Aa^) (^z/(l + /i)"^a + (l + /x)a"f)|0) = 1 + /x. (4.15) 

Consider, for instance, the sum corresponding to F2k,2i- If some term of this sum contains 2j 
pairings of type ()4.15|1 (it is easy to understand that this number should be even and satisfy < j < 
min{/c, /}), then it also contains k — j pairings of type ()4.13|) and / — j pairings of type ()4.14p . The 
total number of such terms is equal to 



Ci^ X C^^j X (2j)! X 



Simplifying this combinatorial factor and taking into account the above remarks, we obtain a general 
formula for the even-even form factor: 

^2^.2'= E io^uT ^ui mf-'{l + ^^?'{^/2V■ (4.16) 

^ (2j)!(A;-j)!(/-j)! 
Analogously, the odd-odd form factor is given by 

min{A;,/} , i M /o; 1M 

^--Mm= E r9- fi^!r;^ ^^r/ ^^ + f^?^^' (^2^ ■ (4.17) 

^ (2j + l)!(A;-j)!(/-j)! 
One may also prove these results by induction, taking as its first step the obvious formulas 

F,,,o = {Ola^'PM = ^ (^)' , i^o,2. = m. [a^f |0) = , (4.18) 

and applying at the next steps the recursion relations 

Fk+i,i = — Fk^i+i + I ( I + fJ, — — — 

i + /i V 1 + /"/ 

Fk,i+i = Fk+i,i + k ( 1 + fx - -^^] Fk-i,i . 

Another problem to be considered in this section is the calculation of correlation functions of 
Bogolyubov fields in the hamiltonian approach. In other words, we want to compute the vacuum 
expectation values of the ordered products of time-dependent operators 

d,,,,,ix)=e-^^dx,,,.e^^, (4.19) 
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where the operator Ox^^^y is defined by (j4.6j) . Let us start with the two-point correlation function 
(C^Ai,Ati,i/i(cti)C^A2,/i2,i'2('^2))- Assuming that a2 > ai and applying the formulas ()2.6p and ()4.18p . one 
obtains 



\2k) (2A;|Ca2,/.2,z.2|0) e" 



2fcm(a2— ai) 



fc=0 



fc=0 



(A;!)2 ^4 ' ' 



fc=i ■ ^ ^ 



The function characterized by such series had already appeared in the computation of the one-point 
function {Oy)- From the comparison of the last expression with the expansion ()4.10p it follows that 



{Ox 



1 - i^iAae-^™^"^-"^) 



-1/2 



(4.20) 



Remark. In order to obtain two-point correlation function of fields, corresponding to (5-interactions, 
one should make in ()4.20p the substitution 

Vi/2m 1/2 /2m 



l + Vi/2m' l + V2/2m' 

and take into account the one-point functions ()4.8p . The result is then given by 



-1/2 



g-2m(a2-ai) 

2m ) \ ' 2m J Am? 
One may easily check that it coincides with the formula (|3.6|) . specialized to the case N = 2. 



(4.21) 



Though form factor series for the multipoint correlation functions of Bogolyubov fields have a 
more complicated structure, compact expressions for these correlators may also be found. Note first 
that the form factors of the operator Ox,^i,u{x) are obtained from the form factors of Ox,^,u by the 
substitution 

Next, the product of any two such operators is again an operator of the form (j4.6p . multiplied by a 
constant. More precisely, one has 

where the parameters A3, /is and z/3 are determined from the comparison of induced linear transfor- 
mations, 

, , , , (1+Aii)^ Aii + Ai2 + /iiAi2 + i^iA2 , (l + /i2)^ 

As = Ai + A2 — , /X3 = 7— TT— ^ , 1^3 = 1^2 + Vi — , 

l-VlM (1 + /ii)(l + /i2) l-Ui\2 

and the coefficient C12 is given by 

C12 = (0|Oai,mi,-i(^A2,^2,.J0) = [1 - Ur\2]-"\ 

Successively applying the above observations, one may reduce any product of time-dependent oper- 
ators to a single operator of type ()4.6|) . multiplied by a constant factor. This constant clearly gives 
the correlation function we are looking for. 
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5 Tau-function of the Schroedinger operator 
with point interactions 

It was shown in Section 3 that the correlation function of certain Bogolyubov fields can be formally 
expressed by the formula ()3.5j) through the (Weinstein-Aronszajn) determinant of the Schroedinger 
operator with 5-interactions. Therefore, it is natural to assume that this correlator has a geometric 
interpretation. In this section, we construct a geometric invariant of the Schroedinger operator, 
simply related with the above correlation function. 

Let us choose a collection a = (ai, . . . , oat) of iV distinct points on the real line and suppose for 
definiteness that ai < 02 < . . . < a^- Consider the A^-parametric family of self-adjoint extensions of 
the operator La, which correspond to separated boundary conditions (that is, to delta-like potentials). 
A general element of this family. Lay, is given by 



rf2 



domLay = {fe H\R) n H\R\a) \ f{a, + 0) - f'{a, - 0) = V,f{a,), (5.1) 

V^eR, j = l,...,N}. 

Self-adjointness of this operator implies the following important statement: 
Proposition 5.1 Ifrn^E C\R, then the Schroedinger equation 

(L,,v + m2)^ = (5.2) 

has no solutions in the domain of Lay. 

Let us now study the spaces of boundary values of certain local solutions of the equation ()5.2j) . It 
will always be assumed that G C\R and Rem > 0. We isolate the points ai, . . . , ajv in the union 
S = UjLi of ^ disjoint open intervals Sj = {Xj,x^), chosen so that Qj G Sj for j = 1, . . . , A^. 
The set S is fixed once and for all, while the coordinates ai,. . . ,a]sf are allowed to vary provided 
each aj stays in Sj. The whole family of operators Lay that satisfy this condition will be denoted 
by Cs- It depends on 2A^ parameters, including the positions of delta-interactions {0^}^^^ 
their strengths {Vj}^^^ j^. Define an auxiliary map 



7T : domL„ v ^ W = ® . . . ® C'^ 



2n times 



where 



Consider now the space of functions from the domain of Lay & Cs, which solve the equation 
()5.2|) in the exterior of S. We will denote by W^^^ the image of this space in W under the map tt. 
Note that the subspace W^^^ (ZW is actually independent of the choice of {aj} and {Vj}. Since for 
any h G W^^^ there exists a function ip G domL^y, such that h = Ti{ip) and {Lay + m'^)ip = on 
R\S, the coordinates of h should satisfy the relations 

ht = C = 0^ (5-3) 
hl^}^ = w,hf^_, h%^=w,h^l'^^\ ^, = e-(-f+i--"), ^ = l,...,Ar-L (5.4) 
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One can define along the same lines the space W^*^ of boundary values of functions from the 
domain of Lay., satisfying the equation ()5.2j) in the interior of S. It is straightforward to check that 
the coordinates of any vector g G W™^ verify 

Jr)=A^«)(|f ). 'V.W) = f f'f^n . « = l.....iV. (5.5) 

where we have introduced the notation 



^ ' l + Vi/2m ' ^ ' l + Vi/2m ' ^ ' 

—mix^—x^) 

m^) = l^m = . ^ , ; ,^ ' . (5.7) 

1 + Vij Zm 



Proposition 5.2 The subspaces W^^^ and W^y are transverse in W . 

■ The subspaces W^^^ and W^y have zero intersection. In the opposite case one would be able to 
correspond to any nontrivial vector w G W^^^P^W^y a global solution of the equation ()5.2|) . However, 
the existence of such solutions is forbidden by the Proposition 15.11 Observing that dim W^^^ = 
dim W^f = 2N, one may now conclude that VT^^* U W^"^^ = W. D 

We can thus associate to any operator Lay & Cs ^ point W^y in the grassmannian Gr{2N,4:N) 
of 2A^-dimensional subspaces of W ~ C^^. Already at the present stage one could define the tau 
function as an invariant of four points of this grassmannian. Such construction, which can be thought 
of as a generalized cross-ratio, has been proposed in jH]: 

Definition 5.3 Let W be a 2k- dimensional complex vector space, and let Gr{k, 2k) denote the grass- 
mannian of its k-dimensional subspaces. Given four points Wi,W2,W3,W4 G Gr{k,2k) in general 
position, the tau function t{Wi, W2, W3, W4) is defined as 

det (Wi ^ W2) 

t{W,, W2, Ws, W,) = ) (- , 

det (Wi ^ W2) 

where Wi — ^ W2 denotes the projection ofWi on W2 along Wj (j = 3,4). 

We will give a more specialized definition of the tau function. First, remark that our grassmannian 
has a distinguished point W^^q, corresponding to Friedrichs extension of the operator La, 

i. e. to ordinary second derivative operator on if^(R). This point can be used to construct the 
det*-bundle over Gr(2A^,4A^). Its fiber at the point U G Gr(2A^,4A^) is a line X*{U) (g) XiWp"^), 
where X{U) denotes maximal exterior power of the space U. 

Another distinguished point of the grassmannian is given by the subspace W'^^'^. The splitting 

W = VT"^* © Wp^ allows to introduce the projection U ^ Wp^ for any U G Gr{2N,AN). This 
projection induces a linear map from X{U) to X{Wp^), and thus defines a canonical section a of the 
det*-bundle: 

a: U^det(u ^ Wp'^ G X*{U) ® X{Wp'). 
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Given another map F : U ^ Wp^*, one would be able to construct in analogous manner a 
trivializing section 6 : U ^ det F, and to define the tau function 

, , det (U !=: Wp') 

It is easy to see that this definition is independent of the choice of bases of U and W^p^. Therefore, 
it indeed gives a function on the grassmannian and, consequently, on the family £5 of Schroedinger 
operators. 

In order to define the map F, let us introduce the notion of auxiliary projections Ps {j = 
1,. . . ,N). Suppose there are no delta-interactions at all and consider a single open interval S' G 
M. The space W{S') = © is decomposed as above into the direct sum of two-dimensional 
subspaces Vr*"*(S") and W'^^^{S'), consisting of boundary values of if^-solutions of Schroedinger 
equation without point interactions on S' and M\S', correspondingly. We will denote by P5/ the 
projection of W{S') on 14^™*(S'') along W^^^{S'). The map F : U ^ is now defined as the 

restriction to [/ of a direct sum of such projections: F = (Pg^ © ... © Psn) 

to check that for any vector ip & U we have 



. It is straightforward 

u 



= Ps,^P^^^ = { j © iV,(0) ( ] , J = 1, . . . , iV. (5.9) 



Let us now obtain the explicit form of the map F and of the canonical projection P : U — > Wp , 
putting U = W^^y. First one should make some choice of coordinate bases. Remark that abitrary 
vectors / G W^^y and g G Wp*^ can be written (as elements of W) in the following way: 

/(i) = /. © N^iVj)f, , g^^^ = ~g, © iV,(0)^, , j = 1, . . . , iV, (5.10) 

where 



h 




and the matrices {NjiVj)}.^^ ^ are defined by the formulas ()5.5|) - ()5.7|) . Hence we can represent / 
and g by the columns 

/=(/7.../^)^ 9={9l...~g],f. 

It follows from ()5.9|) - ()5.10|) that the map F is given in these coordinates by the identity matrix. 

In order to find the representation of P in such coordinates, one should be able to decompose 
any vector / G W^y as f = g + with g G and h G W^^K Let us obtain the relation between 
/ and g. Since the coordinates of h = f — g should satisfy the conditions ()5.H|) - ()5.4|1 . one has 

_ f(l) JN) _ AN) 

yL- ~ J L,-y yR,+ ~ J R,+y 
^ - at = f^:^ -^^ft^ at - & = /gl - fl^^^ J = 1, . . . , AT - 1. 

Using the representations ()5.10|) . we may eliminate from these relations the extra variables /jj'L, Jl^j 
g^^^_, and gp \ {j = 1, . . . , N). Obtained system of linear equations has the form 

{l + Ma,v)f = il + Mo)g, (5.11) 
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where 2N x 2N matrices Ma v and Mq are defined as 



( gi(V2) 

Ti(Fi) Q^iy-i) 

T2(F2) 



V 














Mo 



/ Qi(0) 

ri(o) g2(o) 

T2(0) 



\ 


Qn-M 

T^-i(O) j 



\ 

Auxiliary 2x2 matrices and Tj(Vj), entering these formulas, are given by 











The system (15.111) implies that the tau function ()5.8|) . evaluated at the point W'^^y G Gr{2N, iN), 
is equal to the ratio of determinants 



^ "'^^ det(l + Mo) 



(5.12) 



However, this expression may be simplified. Since Tj{0)Qj{0) = 0, the matrix 1 + Mq can be 
represented as a product of a lower triangular and an upper triangular matrix with identities on 
their diagonals. Therefore, the determinant in the denominator of (|5.12j) is equal to 1. Moreover, 
one can show^ that our tau-function does not depend on the choice of localization, i. e. on the 



coordinates <i Xj'^ \ . In particular, we may put x 



' j=l,-,N 

following representation: 



X 



R 



a J for j = 1, 



N and obtain the 



int \ 

ay) 



det 





1 


Qi 










Ti 


1 


Q2 












1 












Qn 


V 










Tn-i 1 



(5.13) 



where 



-m{a 



'j + l-aj 



a,) 



1 + 



2m 



2m 





-m{ajj^i—aj) 

l + ii 

^ 2m 





2m 



3 



N-1. 



Finally, it is worth mentioning that the tau function T{Wa^y) and the correlation functions of 
Bogolyubov fields, associated to (5-interactions, are related by 



r(w^r^) 



(5.14) 



^The proof is based on Krein's formula for the resolvents of self-adjoint extensions of La and lies beyond the scope 
of this paper. 



14 



This formula may be checked exphcitly for small values of by comparison of (|5.13p and (|3.6|) . 
Example. For the two-point tau function one has 

Comparing this expression with the formulas ()4.8p and ()4.21|) . we may check the validity of ()5.14p 
for N = 2. 



6 Discussion 

It is well-known that in order to describe delta-interactions in two and three dimensions, one needs 
to renormalize the strengths of point sources. This obstacle complicates the construction of fields, 
associated to delta-interactions in higher dimensions. In particular, a naive generalization, Ov{t, x) = 
exp { — I V0'^{t, x)}, does not work, since even the vacuum expectation value of such an operator is 
infinite. It would be interesting to understand how the renormalization can be described (if it indeed 
can) in terms of Bogolyubov fields. 

It would also be interesting to consider instead of the Schroedinger operator the hamiltonian of 
Dirac fermions in 2 + 1 dimensions in the background of magnetic vortices. Its resolvent is the prin- 
cipal ingredient in the computation of induced fermionic vacuum quantum numbers. It also contains 
the information about the scattering of fermions, their bound states, etc. It seems that this resol- 
vent may be expressed along the lines of the present paper (at least, for some values of self-adjoint 
extension parameters) through the ratio of certain correlation functions in the two-dimensional eu- 
clidean quantum field theory, described by the Dirac action with a parity-breaking term. Thus one 
needs to find form factors and correlation functions of the Bogolyubov fields, associated to magnetic 
vortices. The absence of parity-breaking term means that the resolvent is calculated at zero energy. 
Corresponding Bogolyubov fields reduce in this case to the exponential fields of the SG//-theory. 

Yet another area of applications of Bogolyubov fields lies in quantum optics, where they corre- 
spond to the so-called "squeezed" operators jlj, describing the minimum- uncertainty states. 

A Gaussian integration 

Suppose we are given: 

• two real vectors ip = {ipi, . . . , (Pm)'^ and fj, = (/zi, . . . , [inY] 

• a vector J = (Ji, . . . , JmY] 

• two real invertible symmetric square matrices A and B of the corresponding sizes M x M and 
N yiN] 

• an M X N matrix C. 
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If we denote Vip — dipi . . . dipM, T^l^ — dl^i ■ ■ ■ dl^N^ then the following formulas hold: 



Za 
Z 



I 



(27r) 



M+N 



detAdet(5-C^A-iC) 



(A.l) 



M+N , I 

= (27r)^/ Vdet B det {A - CB-^C^) , 
= Af^ + \a-^C {B - C'^A-^cy^C'^A-^ 



(A.2) 



Now let the vectors </? and be complex, and let the matrices A and B be hermitian. Then, 
introducing the measures 



M 



N 



VipVCp = 2^ d{Re ipj) d{lm ipj), 

one may obtain the analogs of the above relations: 

(27r)^ 
det^l ' 



Za^ j VipVip e"'^^^'^ 
Z^^ I VipVCpe 



Z-/ j V<pV(p ^.^^e-^'^^^A'l 



Zj^ \ Vip Vip (pi(pj e 



-ifi'^Aip+j'^'P+ifi^ J 



Za,b,c = j Vip Vip V^Vfi e-^'^^-^'^^-^''^^-^''^'^ = (27r)"'+^/det 5 ) = 

= (27r)^+^/[det A det {B - CU-^C)] = (27r)^+^/[detS det {A - CB-'C^)] , (A.3) 



Za~b,c J ^/^ e-'^^^'^-'^^-^'^-'^^^'^-'^^^^'^ = 

{A - CB-^C^)~I = A-^ + \a-^C {B - C^A-^C)~^C^A-^ 



(A.4) 
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